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ABSTRACT 


The Nysa-Polana complex is a group of low-inclination asteroid families in 
the inner main belt, bounded in semimajor axis by the Mars-crossing region 
and the Jupiter 3:1 mean-motion resonance. This group is important as the 
most likely source region for the target of the OSIRIS-REx mission, (101955) 
Bennu; however, family membership in the region is complicated by the presence 
of several dynamically overlapping families with a range of surface reflectance 
properties. 

The large S-type structure in the region appears to be associated with the 
parent body (135) Hertha, and displays an (ep,ap) correlation consistent with a 
collision event near true anomaly of ~ 180° with ejecta velocity Uej ~ 285 m/s. 
The ejecta distribution from a collision with these orbital properties is predicted 
to have a maximum semimajor axis dispersion of Saej = 0.005 ±0.008 AU, which 
constitutes only a small fraction (7%) of the observed semimajor axis dispersion, 
the rest of which is attributed to the Yarkovsky effect. The age of the family is 
inferred from the Yarkovsky dispersion to be OOO^gQ My. 

Objects in a smaller cluster that overlaps the large Hertha family in proper 
orbital element space have reflectance properties more consistent with the X-type 
(135) Hertha than the surrounding S-type family. These objects form a distinct 
Yarkovsky “V” signature in (op, H) space, consistent with a more recent collision, 
which appears to also be dynamically connected to (135) Hertha. Production of 
two families with different reflectance properties from a single parent could result 
from the partial differentiation of the parent, shock darkening effects, or other 
causes. 

The Nysa-Polana complex also contains a low-albedo family associated with 


(142) Polana (called “New Polana” by Walsh et ah 2013), and two other low- 


albedo families associated with (495) Eulalia. The second Eulalia family may be 
a high-ap, low-ep, low-fp component of the hrst Eulalia family-forming collision, 
possibly explained by an anisotropic ejection held. 


Keywords: Asteroids; Asteroids, dynamics 
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1. Introduction 


The low-inclination, inner main belt asteroid population is an important source region 
for near-Earth objects (NEOs) and terrestrial planet impactors in the past, present and future 
(Scholl and Froeschle 1991 Morbidelli and Vokrouhlicky 2003). This region, which includes 
the Nysa-Polana complex, poses a challenge for family identihcation and analysis techniques, 
because it is dynamically very crowded, with many families overlapping in proper orbital 
element space. Recent efforts have used newly available color information from the Sloan 
Digital Sky Survey’s Moving Objects Catalog (SDSSMOCpI [Ivezic et al.||2002 ) and albedo 


information from the Wide-held Infrared Survey Explorer (WISE/NEOWISEQ 


Wright et ah 


2010 Mainzer et al.|20lT) to augment the orbital databases in an effort to further distinguish 


families on the basis of their unique rehectance properties (|Parker et ah 2008 Masiero et ah 


2013 Carruba et al. 2013 Milan! et ah 2014; Walsh et ah 2013 Dykhuis et al. 2014). 


Analysis of the taxonomically diverse and dynamically overlapping families within the 
Nysa-Polana complex has suffered from a long and complicated history of identihcation and 
re-identihcation, often with changes in membership, taxonomy, and nomenclature. In order 
to place the current work in context, we hrst review the history of these identihcations 
and the rationales that led to each. The various names assigned to each family structure 
throughout the literature are listed in Table In the interest of connecting our work to 
future ehorts, we also identify the Nysa-Polana region with the Family Identihcation Number 
(FIN) of 405, assigned to the region by Nesvorny et al. and Masiero et al., both submitted. 
In this classihcation scheme, the Herthal family would correspond to 405a, the Polana family 
to 405b, and the Eulalial family to 405c. The IIertha2 and Eulalia2 families would await 
further conhrmation from additional data. 


The Nysa-Polana complex was hrst identihed as a single family by Brouwer (1951). 


Zellner et al. (1977) subsequently distinguished two family structures in the region, and 


suggested an origin scenario that attempted to include both of the largest neighboring objects 
(44) Nysa and (135) Hertha as fragments from the same parent body. Williams (1979) 
conhrmed the existence of two families, labeled as ‘W-24” and ‘W-160.” With the addition 


of rehectance property data from the Eight-Color Asteroid Survey (EGAS), Tedesco et al. 
(1982) discovered several objects in the so-called “Nysa” family to be the otherwise rare 
“F” spectral type (using the taxonomy scheme developed by the authors for the EGAS). 
This spectral dissimilarity led Bell (1989) to conclude that the “Nysa” family members were 


Tttp: //www.astro.washington.edu / users / ivezic / sdssmoc / 
^http://irsa.ipac.caltech.edu/Missions/wise.html 
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unrelated to the E-type (44) Nysa, and should rather be associated with the F-type (142) 
Polana, labeling (44) Nysa an interloper in that family. Kelley et ah (1994) suggested an 


origin scenario involving partial differentiation in attempt to link the E-type (44) Nysa and 
M-type (135) Hertha as “parents” of the F-type family, based on a shared silicate absorption 
feature near 0.9 pm. 


This work 
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Abbreviations: Ze77 = 

= jZellner 

etallllTg^; Wi79 

= [William^ B89 

= |Bell|(|1989|l; Za95 =|Zappala 

et al.|l|1995 


COl = |Cellino et al.| ( |2001[ |; MD05 = |Mothe-Diniz et al.|p005^ ; D12 = |Delbo’ et al.| ( |2012^ ; Wal3 = [Walsh et al.| ( |2013^ ; 
M14 = |Milani et al.|(2014| |. 


Table 1:: Summary of the various names of Nysa-Polana complex families used in the literature, showing 
sources that modihed the family nomenclature. See the text for a more detailed summary of each literature 
source. The first column shows the five structures of the Nysa-Polana complex, using the naming conventions 
of the present work. The spectral types in column 2 are broad categories only, meant to reflect the colors 
observed by the SDSS; the individual literature sources often use subcategories or parallel categories (such as 
Xc, Xk, E, F, B, etc.). Ditto marks indicate that the Eulalial structure, defined in Section 4.2 was typically 
joined to the Polana structure until its distinction in the work of Wal3. The question mark expresses low 
confidence reported by Wal3. MD05 lumped the low-ep portion of our Herthal structure together with our 
Hertha2 structure, retaining the name Mildred for the higher-ep portion of Herthal. 


Zappala et al. (1995) used two different clustering techniques to identify family members 


in the Nysa-Polana region, and found two overlapping families that the clustering methods 
failed to reliably distinguish. The Polana family was labeled as a possible subsequent collision 
on a Nysa family member. 


Cellino et al. (2001) undertook a spectroscopic campaign in an attempt to clarify mem¬ 


bership in the region. They found three objects with X-type taxonomy among the F-type 
Polana objects and S-type Hertha objects. They renamed the Hertha objects as the “Mil¬ 
dred” family, based on the dissimilarity between the M-type (135) Hertha and the S-type 
family previously associated with it. 


Further observations and dynamical analysis by Mothe-Diniz et al. (|2005) and Alvarez 


Candal et al. 

(2006 

) identihed three more objects as X-types, including (44) Nysa and (135) 

Hertha. 

Mothe-Diniz et al. 

(2005 

) separated the Mildred family into two groups based on 


their eccentricities, and noted that the lower-ep group contained a mix of X- and S-type 
objects, while the higher-ep group contained primarily S-type objects only. They combined 
the lower-ep portion of the Mildred family with what we now call the Hertha2 family to 
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produce a “Hertha” family, while the higher-ep component retained the name “Mildred.” 
In addition, they renamed the Polana family as “McCuskey.” More recently, |Delbo’ et ah 


(2012) tested a new classification algorithm on the region, and identified three groups of 


objects based on reflectance properties, confirming the presence of S-, X- and C-type spectra 
in the Nysa-Polana region. 


Walsh et al. (2013) used the recent data available from WISE/NEOWISE to study 


the dark (B- and C-type) objects in the Nysa-Polana complex, which until then had been 
considered a single family (usually called “Polana”). They identihed the structures of at least 
two old families, which they associated with (142) Polana and (495) Eulalia. A third structure 
at lower eccentricities was identihed with (112) Iphigenia, although with low conhdence in 
its signihcance or its independence from the Eulalia family. 


The most recent analyses (Masiero et al. 2013 Carruba et al. 2013; Milan! et al. 2014) 


incorporated major results from large rehectance surveys, primarily WISE and SDSS, and 
used automated clustering methods applied to this extended parameter space. In general, 
they conhrmed the presence of at least one low-albedo (B- or C-type) collisional family and at 


least one high-albedo (S-type) collisional family in this region. Milan! et al. (2014) identihed 


the hrst cluster in Table with a diherent parent asteroid, (3583) Burdett. 

Here, in addition to seeking clusters in the dynamical and rehectance parameter space, 
we seek signatures of dynamical evolution that distinguish multiple families within a sin¬ 
gle cluster. Specihcally, by considering cluster distributions in ap — Ef space, the distinct 
“V-shape” patterns of Yarkovsky evolution can be recognized. This technique was demon¬ 
strated recently for the “Eulalia” and “New Polana” families by Walsh et al. (2013), and was 


also successfully applied in our earlier analysis of the large and dihuse Flora family, which 
dominates the inner main belt ( Dykhuis et al.||2014 ). In addition, for the Herthal family, we 
extract collisional parameter information from the slope of the family’s ap-ep distribution. 

In Section we give an overview of the geography of the Nysa-Polana region in proper 
elements and rehectance properties. Then, we separate the objects on the basis of their a* 


color parameter (as dehned by Ivezic et al. 2001), and study hrst the objects with a* > —0.015 
(Section]^, which we show to be one of two families likely associated with (135) Hertha 
(designated the Herthal family, see Table [^. Next we study the objects with a* < —0.015 
(Section]^, identifying four structures and their likely parent objects. 
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2. Overview of the region 


The orbital distribution of asteroids in the Nysa-Polana region (2.1 AU < ap < 2.5 AU, 
0.11 < ep < 0.22, and 0.02 < sinzp < 0.07) with reliable proper orbital element data is 
shown in Fig. The catalog of synthetic proper elements is the same as that used in |Dykhuis 
|et al.| ( [2M^ , computed using the OrbitO software available from the Asteroids Dynamic Site 
(AstDyS)^ The blue and red boxes in Fig. la roughly represent the regions associated with 
the low-albedo (C-type) and higher-albedo (S-type) groups, respectively, that have typically 
been identihed in the region. Asteroid (142) Polana lies within the former range, (135) 
Hertha in the latter, and (44) Nysa in neither. 


EDITOR: PLACE FIGURE □ HERE 


Figs, lb and Ic show the distributions in proper semimajor axis, ap, vs. absolute 
magnitude, H, of the objects from the sub-regions marked in Fig. A “V” shape in 
ap — H space is generally attributed to the Yarkovsky effect, which yields a size-dependent 
spread in semimajor axes for a collisional family. The signatures of several overlapping V’s 
are evident in Figs. [TO and suggesting the presence of multiple families. 


The boundary of each V-shaped Yarkovsky envelope can be described by the limiting 
distance in semimajor axis Aap from the family’s central semimajor axis location, which is 
generally assumed to be the proper semimajor axis of the largest remnant of the family. The 


distance is described by the relation (from Vokrouhlicky et ah 2006): 


Aap = C- 10^/^ (1) 

The O parameter (which for convenience we later refer to in units of mAU = 10“^ AU) is 
directly related to the family’s dynamical evolution due to the Yarkovsky effect. Families 
with an older age will show more semimajor axis evolution, and reach higher boundary values 
of C than their younger counterparts; thus the boundary of the C parameter distribution 
for a given family can be used as a simple means to estimate the age, assuming appropriate 
calibration of the Yarkovsky semimajor axis drift rates and estimation of the semimajor axis 
dispersion due to the original collision ejection field. 

The largest objects in the Nysa-Polana region — (44) Nysa, (135) Hertha, and (142) 
Polana — all lie quite close in semimajor axis to the Mars 1:2 mean-motion resonance (dotted 


^http: //hamilton.dm.unipi.it/astdys/ 
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line in Figs. 

lb 

and 

Ic 

, located at 

Op = 2.419 AU). The effects of this resonance were studied 

in detail by 

Gallardo et ah 

(2011 

), and must be taken into account in studies of the families 


in the inner main belt, as discussed in Section]^ 


Color information for the 8396 Nysa-Polana region objects that have SDSS color data 
(but not necessarily albedo data) are plotted in Fig. Two clusters of objects are apparent 
in this plot, roughly associated with the S-type (high a* color) and C-type (low a* color) 
spectral categories. In Section we consider the cluster with a* > —0.015, which is a 
recognizable dynamical family. The objects with a* < —0.015 represent multiple dynamical 
families, as shown in Section This multiplicity is already indicated by the presence of 
a number of X-type objects in the region with a* colors intermediate between the S- and 
C-type ranges; in Fig. [^they appear as the subtle grouping between —0.1 < a* < —0.0. 

EDITOR: PLACE FIGURE [2] HERE 


3. High-a* objects: The Herthal family 
3.1. Family identification and characterization 

The proper element distribution for the Nysa-Polana region objects from Fig. with 
a* > —0.015 is shown in Figs. and These hgures suggest that the high-a* structure 
(within the red box in Fig. is associated dynamically with parent object (135) Hertha. 
We designate this cluster as the “Herthal” family (cf. Table for previous designations in 
the literature), with the “1” added because (135) Hertha is also associated with a smaller 
cluster of objects (“Hertha2,” described in Section]^. 

EDITOR: PLACE FIGURE [3] HERE 


The distribution within the red box in Fig. could be interpreted as two separate 
clumps around ep, sin ip of (0.170, 0.045) and (0.185, 0.040). These two clumps were named 
the “Hertha” and “Mildred” families, respectively, by Mothe-Diniz et al.| ( |2003 ), cf. Table 
However, the op-cp distribution of the objects within the red box (Fig.|^ shows a continuous 
and uniform trend of lower ep with increasing ap, consistent with the ejection held from 
a single large-scale collision event. Moreover, as we show below, the Yarkovsky V shape 
is similarly continuous within and between the two populations. Thus we interpret the 
structure within the red box in Fig. |^as a single collisional family, Herthal. 


EDITOR: PLAGE FIGURE i HERE 
















The absolute magnitude distribution of the objects of the Herthal family yields informa¬ 
tion about its dynamical evolution under the influence of the Yarkovsky effect. The ap — H 
distribution plotted in Fig. shows the V shape characteristic of a collisional family evolved 
under the Yarkovsky effect. The rightmost edge of the V shape is truncated by the Jupiter 
3:1 mean motion resonance (represented by the blue dotted line), and the objects at high 
eccentricities and low semimajor axes are affected by encounters with Mars, modifying the 
distribution from an “ideal” V shape. The low-op, low-ep objects evident in Fig. |^are here 
shown to be primarily interlopers in the family, included in this sample only on the basis of 
their inclusion within the family ranges in ep and sin ip and a*. 

EDITOR: PLACE FIGURE [5] HERE 

In principle, the age of the family can be interpreted from the distribution of C pa¬ 
rameters of the family members (Equation]^, the outer boundary of which describes the 
V-shaped envelope in ap vs. H. However, the correlation between ep and ap among the 
Herthal family objects (Fig. means that the center of the Yarkovsky V shifts to lower 
ap for larger eccentricities. This complicates the identihcation of the Yarkovsky envelope 
that best hts the boundary of the family in semimajor axis; none of the example Yarkovsky 
envelopes plotted in Fig. [^match the boundary well. The full semimajor axis width in Fig. 
is due to the skewed ep vs. ap distribution in Fig. which is related to the initial collision 
rather than the subsequent Yarkovsky drift. 

This skew must be removed prior to interpretation of the family’s age via the observed 
Yarkovsky spread in ap. To do this, we need to adjust each member’s semimajor axis by an 
amount 


ap — ap — — ■ (ep — epo), (2) 

m 

where m represents the slope of the ap vs. ep distribution, and ep^o is the family center in 
eccentricity. Given appropriate values for m and ep^o, we can apply a shear transformation 
to the distribution of objects to remove the eccentricity dependence and determine the 
Yarkovsky-induced spread in ap. 

In order to hnd the slope m of the observed distribution, we must take into account 
the effects on the shape of the distribution due to the various resonances in and around 
the Nysa-Polana region (Fig. |^. The Mars-crossing region removes low-ap, high-ep Herthal 
objects; the Jupiter 3:1 mean-motion resonance removes high-ap, low-ep objects; lastly, the 
Mars 1:2 mean-motion resonance depletes family members both below ap = 2.43 AU and 
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above ep = 0.195 as they drift outward in semimajor axes due to the Yarkovsky effect and 
exit the region before drifting back. 

The only “pristine” edge of the family that remains untouched by resonances (and 
thus still reflects the semimajor axis dependence of ep from the original collision) is the 
lower left edge, at 0.15 < ep < 0.185, where the sloped boundary of the family is visible 
(Fig-i- We use this boundary to estimate the slope of the family; specihcally, we hnd the 
boundary at each value of eccentricity within the range 0.15 < ep < 0.185, using a stepsize 
of 6e = 0.001 and an eccentricity width of Ae = 0.01. The latter parameter specihes the 
window in eccentricity which dehnes the sample of objects used to hnd the boundary at 
each eccentricity step. The boundary is dehned for this purpose as the point of maximum 
increase in density of objects, or a maximum in the derivative of the density function (see 


Appendix A.l of Dykhuis et ah (2014) for discussion of the kernel density estimation and its 
derivatives). A linear £t to the boundary values yields a slope of m = (—0.50 ± 0.04) AU“^, 
represented by the dotted line in Fig. 

The slope m is a signature of the semimajor axis and eccentricity dispersion of the 
Herthal family that resulted from the original collision. In particular, m is determined by 
the true anomaly, /, of the parent asteroid at the time of its disruption. This relationship 
is described by Gauss’ equations, which give the change in orbital elements for a collision 
fragment that experiences a velocity impulse An (see, e.g. Zappala et al.||2002 ). Gauss’ equa¬ 
tions describe the behavior of the collision fragments in unperturbed (osculating) elements; 
in the Hertha region, the proper elements show similar behavior. 

We plot the proper element distributions for several different values of / in Fig. 


(cf. Nesvorny et al. 2002, Fig. 1), along with their corresponding slopes, m. Each ellipse 
represents an isotropic ejection held with An = 285 m/s, with An chosen to approximate 
the eccentricity range of the family, estimated by the noticeable decrease in number density 
near ep = 0.15 and ep = 0.205. The ellipses were generated from 100 test particles, with 
osculating orbits calculated via the Gauss equations, and proper elements calculated via two- 
million-year integrations using OrbitO. The ehects of the Mars 1:2 mean-motion resonance 
at Op = 2.419 AU are observed in the proper element calculations. The lowest possible 
value for m occurs at / = 180°, where m = —0.50 AU“^. Thus, assuming the original 
ejection velocity held was fairly isotropic, it appears the collision that formed this family 
occurred when the parent body was near aphelion, with a true anomaly within the range 
170° < / < 190°, calculated using the uncertainty on the slope as determined from the 
boundary values above. 


EDITOR: PLAGE FIGURE [6] HERE 
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In order to apply the correction (Eq. to remove the effect of the initial slope in the 
(ep,ap) distribution, we need to evaluate cp^o in addition to m. The value ep^o is the proper 
eccentricity at the family’s center, which we take to be the barycenter of the population. To 
minimize the influences of the nearby resonances on the barycenter, we consider only the 
largest objects (12.0 mag < H < 14.0 mag), which have experienced the least amount of 
Yarkovsky semimajor axis dispersion toward significant resonances. We further restrict the 
sample to 2.31 AU < op < 2.5 AU, 0.136 < ep < 0.23, and 0.033 < sin ip < 0.051, in order 
to conservatively include all Herthal family members while largely eliminating members of 
the nearby Massalia and Flora families. We also removed from the sample all objects with 
albedos lower than 0.14, which are likely to be associated with the overlapping darker families 
discussed in Section This does not remove all of the contamination from those families, 
due to the incompleteness of the WISE database; however, such contamination would likely 
be primarily from smaller objects, and thus would have little effect on the determination of 
the barycenter of the Herthal family. Lastly, we removed asteroid (135) Hertha itself from 
the barycenter consideration, due to the uncertainty surrounding its inclusion in the family. 

The barycenter of the sample described above is op^o = 2.412 ± 0.015 AU and ep^o = 
0.177 ± 0.003, with the uncertainty dominated by the exclusion of (135) Hertha from the 
sample. The corrected ep vs. a'p distribution, with Op calculated using Eq. with m = 
—0.50 AU“^ and ep^o = 0.177, is shown in Fig. 

EDITOR: PLACE FIGURE [7] HERE 


3.2. Post-collision history of the family 

The corrected semimajor axis dispersion of the Herthal family records information about 
the family’s post-collision evolution under the Yarkovsky effect. In general, the most useful 
Yarkovsky drift information is preserved in the outer edges of a family, specifically the extent 
of drift at each size range, as parameterized by C. The age of the family is obtained from the 
value of C for which Eq. gives curves in a'p — H space that best fit the outer edges of the 
distribution. However, the Jupiter 3:1 resonance at high ap and the Mars-crossing region at 
low Op have significantly sculpted the boundaries of the Herthal family (Fig. [^. The former 
renders the high-ap (right) edge of the family unreliable; the Mars-crossing region affects 
the low-ap (left) edge for ep > 0.185. Accordingly, in order to determine an appropriate 
boundary C value, we consider only the objects with ep < 0.185, and fit only the lower-ap 
edge of their distribution in a'p — H space. 

To determine the value of C that best fits that edge, we calculate for each asteroid the 
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value of C that corresponds to its a'p and H values (Eq. [^. Fig. shows the distribution 
of these C values, with the curves representing the kernel density estimate (KDE) for the 
distribution along with its derivatives (see Appendix A.l of Dykhuis et ah (2014) for an 


explanation of the KDE and its derivatives). The best-£t edge of the family is the C value 
where there is a maximum in the second derivative of the KDE (the blue curve in Fig. [^. 


EDITOR: PLACE FIGURE i HERE 


There are comparable maxima in the second derivative at C = —0.045 mAU and 
—0.037 mAU; however, inspection of the distribution in a'p — H space (Fig. shows that 
the latter boundary would lie well inside the dense population, not at its edge. By contrast, 
C = —0.045 mAU (black curve in Fig. does indeed £t the outer edge of the distribution 
quite well. The uncertainty in the C boundary (±0.005 mAU, represented by the dashed 
curves in Fig. is dominated by the uncertainty in the ap^ center of the family. The C 
boundary dehning the edge of the family is much tighter for the corrected a'p — if distribution 
than it was for the uncorrected ap—H distribution plotted earlier in Fig.[^(C' = ±0.045 mAU 
vs. C = ±0.075 mAU). 


EDITOR: PLACE FIGURE [9] HERE 


The value of G = 0.045 mAU fits the leftmost edge of the distribution quite well for 
the larger objects {H < 17 mag), but fails for the smallest objects {H > 17 mag), which 
is not surprising for a couple of reasons. First, the observational bias in the database af¬ 


fects the completeness of the family sample for objects fainter than magnitude 16.5 (Bottke 
et ah 2014). Second, the effects of “variable” or “stochastic” YORP on the drifting family 


members can cause their spin rates to random walk toward the typical YORP end states 
(mass shedding or non-principal axis rotation), which applies a size-dependent “brake” to 
the smallest objects ( Bottke et al.||20f3 , |2014 ). Stochastic YORP simulations of the Eulalia, 
“New Polana” and Erigone families by these authors demonstrated behavior similar to what 
is observed here. 

The calculation of the age of the family from the observed Yarkovsky semimajor axis 
dispersion requires an estimate of the initial semimajor axis spread due to the original colli¬ 
sion itself. This estimate is typically difhcult to obtain, but in the case of the Herthal family, 
we can place an upper limit on the collision spread from the ap — ep distribution. As shown 
in Section 3.1, the slope of m = —0.50 AU“^ indicates a collision near aphelion, which would 


result in minimal initial semimajor axis dispersion in a'p. The maximum expected value of 
dispersion is da^j = 0.005 ± 0.008 AU, which represents the width of the proper semimajor 
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axis dispersion for a collision at true anomaly / = 180°, with uncertainty derived from the 
uncertainty in / (which results from the uncertainty in the slope m). 

The semimajor axis dispersion due to the original collision is thus expected to be a 
small fraction of the Yarkovsky dispersion, resulting in an overestimation of the boundary 
C parameter of less than 0.003 mAU, which is less than the uncertainty on the boundary C 
parameter (0.005 mAU). Hence the initial semimajor axis dispersion has little effect on our 
calculation of the age from the Yarkovsky drift. 


We use the inner boundary of the C parameter distribution to determine the age of the 
family via the relation (Dykhuis et ah 2014): 


t = 


1329 ■ \C\ ■ ap,o^ 
cyA/W(l -Pv) cose 


( 3 ) 


where t is the time (in My) since the collision, ap^o is the time-averaged proper semimajor axis 
of the family members (in AU), e is the obliquity, and C is given in AU. The parameter cy 
contains information about the asteroid’s material properties (thermal conductivity, specihc 


heat, material density), which are not known a priori. As in Dykhuis et al. (2014), we adopt 


the value for cy obtained from measurements of the precession rates of the Karin family, 
which yield a cy parameter in the range of 0.0026 AU^km My“^ < cy < 0.0035 AU^km My“^. 

With values oi C = ±0.045 mAU, albedo pv = 0.25, ap,o = 2.412 AU, and cy = 
0.00305 AU^km My“^, we find an age for the Herthal family of OOOj^gQ My. The error in this 
estimate is dominated by the uncertainty in cy. 


The age estimate also does not account for possible systematic errors due to the assump¬ 
tion of similar material properties between the members of the Herthal and Karin families. 
In particular, differences in bulk densities would affect the cy parameter appropriate to both 
families. The average bulk densities for the two families are difficult to estimate, due to the 
lack of available density data for the family members. The only Hertha family object with 


density data is (135) Hertha itself (Carry 2012); its taxonomic type of makes it a poor 
standard for comparison with the S-type Karin family. The Karin family itself has members 
with density measurements, but the Koronis family, of which Karin is a subfamily, has den¬ 
sity measurements for two of its members, (243) Ida and (720) Bohlinia; their densities are 


2.35 ± 0.29 and 2.74 ± 0.56 g/cvsr, respectively (Carry 2012). These are both consistent 


with the averages for the S-type class, and thus are assumed here to be consistent with the 
unknown Hertha family average. 
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4. Low-a* objects: The Hertha2, Polana and Eulalia families 


We now consider the objects in the Nysa-Polana region with a* < —0.015. The (ep, 


sinzp) distribution of these objects is plotted in Fig. 10 Several structures appear within 


the blue box; the diffuse grouping near (142) Polana and the dense grouping near ep = 0.145 
and sin ip = 0.047 are here referred to as the Polana and Eulalial families, respectively 
(Table [^. The combined families had long been recognized as a single “Polana” family until 
their distinction in the work of Walsh et ah (2013|). 


EDITOR: PLACE FIGURE M HERE 


Comparison of Fig. with Figs. and shows that the Herthal family, roughly 
bounded by the red box, has been removed. However, a smaller, dense cluster of low-a* 
objects lies near the location of asteroid (135) Hertha, within the range of the now-removed 
Herthal family. We refer to this cluster as Hertha2. 


4.1. The Hertha2 family 


The Hertha2 cluster lies very close to asteroid (135) Hertha in ap (Fig. [lT|), as well as in 
ep and sin ip, suggesting a genetic relationship. In fact, the barycenter of the objects within 
the cyan box in Fig. 10 (those plotted in Fig. is at ap = 2.426 AU, nearly equal to the 
value for (135) Hertha itself, ap = 2.429 AU. The cluster displays the characteristic V-shape 
in ap — H indicative of Yarkovsky evolution, further identifying the Hertha2 structure as 
a collisional family. The V-shape envelope corresponds to C = ±0.012 mAU (signihcantly 
narrower than that observed for the Herthal family, C = 0.045 mAU, Section [^. This 
tighter V could represent either a younger age or a higher average density for the Hertha2 
family; the implications are discussed further in Section 


EDITOR: PLACE FIGURE [U] HERE 


The Hertha2 family also appears as a distinct cluster in reflectance properties. Fig. 12 
shows the color distribution similar to Fig. but restricted to the objects within the same 
dynamical range as the Hertha2 family (2.1 AU < ap < 2.5 AU, 0.1690 < ep < 0.1812, and 
0.0420 < sin ip < 0.0471). The Hertha2 family forms a diffuse grouping centered around 
a* = —0.07; the members of the large Herthal family that reside in the same dynamical 
space are distinct from the Hertha2 members, clustering instead near a* = 0.13. 


EDITOR: PLACE FIGURE M HERE 
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4.2. Other dark families 


Other family structures among the dark asteroids in the Nysa-Polana region can be 
teased apart by considering the ap — H distribution in the context of the Yarkovsky evolution 
model, an approach utilized by Walsh et ah (2013). The ap — H distribution of the low- 

This hgure includes for context 


objects in the Nysa-Polana region is shown in Fig. 13 


the cyan boundary C parameter lines for the Hertha2 family (Fig. [IT| ), and plots new C 
parameter lines for the two distinct families identihed by Walsh et ah ( |2013 ) (cf. Fig. 2 of 
that work), which are centered on (142) Polana and (495) Eulalia. These two families are 
visible in Fig. 13; we refer to them as the Polana and Eulalial families, respectively. The 
latter is one of two structures that we identify with (495) Eulalia (the second is discussed 
below). 


EDITOR: PLACE FIGURE M HERE 


Fig. 10 demonstrates a phenomenon that was noticed hrst by Walsh et ah (2013): (495) 


Eulalia is signihcantly offset from its nominal family in ep (ep = 0.12 vs. the family average 


of ep = 0.145). Walsh et ah (2013) modeled the long-term evolution of (495) Eulalia, 
and demonstrated that proximity to the Jupiter 3:1 resonance results in variability of the 
asteroid’s proper eccentricity between 0.11 and 0.15 over 500 Myr. 

In order to explore whether any other asteroids might be related to (495) Eulalia, we 
modify the range under consideration to 0.05 < ep < 0.2 and 0.0 < sinfy < 0.08, and replot 
the distribution of low-a* objects in this new range (Fig. 14). A diffuse grouping appears 
near ep = 0.11 and sinfy = 0.035, indicated roughly by the magenta box in Fig. This 


grouping was noted by Walsh et al. (2013) after removal of the dense Eulalial cluster, and 


was tentatively associated with (112) Iphigenia (indicated by the green diamond). 


EDITOR: PLACE FIGURE M HERE 


Fig. [15] shows the ap — H distribution of the objects in the green and magenta boxes 


in Fig. lA, that is, the objects of the Eulalial and the putative Iphigenia cluster. The latter 
objects show similar structure to the Eulalial objects, and can be described by the same 
Yarkovsky C envelope. A preliminary best fit C envelope for each is at C = ±0.075 mAU 


(slightly different from that found by Walsh et al. (2013), C = ±0.092 mAU). Note that the 
two groups display a distinct offset in semimajor axes (Jap ~ 0.035 AU); this offset, coupled 
with the identical boundary C envelopes for the two clusters, leads us to designate the second 
cluster as “Eulalia2” and interpret it as possible evidence for an anisotropic ejection field 
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resulting from the Eulalial collision event, rather than associate it with (112) Iphigenia. The 
implications of the structures of these two clusters are discussed further in Section 5.2[ 


EDITOR: PLACE FIGURE M HERE 


The semimajor axis offset between these two populations was not discussed by [Walsh 


et ah (2013), and indeed is not apparent in the corresponding hgure in that work (their 
Fig. 20). Walsh et ah relied on the hierarchical clustering method (HCM) to identify 
groupings of objects for comparison. While the HCM is ideal for identifying clusters in 
multidimensional space, it is less ideal for comparison of neighboring families, as it tends 
to include in each sample objects from the neighboring family or families, which muddies 
the comparison. Instead, we select for comparison two distinct sets of objects (within the 


green and magenta boxes in Fig. 14) for which we expect minimal contamination between 


sets, which enables us to identify the small offset in op among the Eulalial and Eulaha2 
structures. 


5. Discussion 


We have identihed hve distinct clusters of asteroids within what has generally been 
called the Nysa-Polana orbital region of the asteroid belt. Two of these clusters appear to 
be associated with collisions on (135) Hertha; two other clusters appear to be associated 
with (495) Eulalia; one is associated with (142) Polana. None of these is associated with 
(44) Nysa. 

This rehned interpretation of the dynamical structure of the Nysa-Polana region has 
been made possible by the availability of color and albedo information from the SDSS and 
WISE projects, which have increased the dimensionality of data available for cluster analyses 
beyond simply the proper orbital element distribution. The analysis of the clusters within 
the Nysa-Polana region depends particularly upon this increased dimensionality, due to the 
overlap within the region of multiple large families with varied reflectance properties. 


Other recent work has explored clusters within this newly extended parameter space 


using an automated approach based on the hierarchical clustering method (Masiero et ah 


2013 Milani et ah 

2014; 

Carruba et ah 

2013 

). Our approach differs from theirs in that we 

(like Walsh et ah 

(2013 

)) have also incorporated the effects of Yarkovsky semimajor axis drift. 


which results in recognizable patterns in the relationship between semimajor axis spread and 
asteroid size (equivalently, the absolute magnitude), specihcally the characteristic V shape 
in Op vs. H. In addition, we use the slope of the op — ep correlation among the Herthal 
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family members as a constraint on the orbit confignration of the parent object at the time 
of the collision, enabling an improved evalnation of the Yarkovsky semimajor axis dispersion 
of the family independent of the collisional dispersion. 


5.1. Herthal and Hertha2 


The collision that formed the Herthal family SOOlgQ Ma (Section distribnted ejecta 
across nearly the entire span of the inner main belt in semimajor axis, depositing material 
directly into Mars-crossing orbits at the time of the event. Under the inflnence of the 
Yarkovsky effect, the Herthal family has since provided a steady stream of material of 
varions sizes into both the Mars-crossing region and the Jnpiter 3:1 resonance. Hence, it is 
plansible that the Herthal family has been a sonrce popnlation for significant impact events 
of varions magnitndes on Earth in the past 300 million years, a time span which inclndes 


the event hypothesized to be connected with the Permian-Triassic extinction (Ranp and 


Sepkoski 1982 Jin et ah 2000, 251.4 Ma,) 


The weak Mars 1:2 resonance affects the semimajor axis and eccentricity distribntion 
of the Herthal family (Gallardo et ah 2011). The effect of this resonance is most notice¬ 


able in the dearth of objects in a depleted zone with op < 2.43 AU and ep > 0.195 in 
Fig-i compared for example with the density of objects jnst below (ep < 0.19) in the same 
fignre. When objects drifting from the family center toward higher semimajor axes reach 
the resonance, it diffnses their eccentricities, resnlting in higher average eccentricity of the 
resonant popnlation. Thns, for objects starting at low eccentricities, the boost moves them 
to a higher-ep section of the Herthal family; for objects starting at high eccentricities, the 
boost may remove them from the family altogether, hence the pancity of bine dots in the 
depleted zone. This pancity demonstrates also the effects of the combined Yarkovsky and 
YORP effects on the Herthal family; objects that wonld otherwise have drifted back and 
forth as a resnlt of YORP cycles, continnally re-filling the depleted zone, have instead been 
removed from the family to higher eccentricities. 

The presnmed parentage and nomenclatnre of the Herthal family has varied thronghont 
the literatnre (Table [^. Unlike the small, low-a* Hertha2 family, whose objects show re¬ 
flectance properties consistent with X-type spectra and thns similar to (135) Hertha itself, the 
large Herthal family shows average reflectance properties consistent with S-type spectra. S- 
type objects wonld not generally be expected to originate from an X-type parent body, given 
the standard pictnre of S-type objects originating from chondritic, nndifferentiated parent 


objects. However, Weiss and Elkins-Tanton (2013) snggested that some chondritic bodies 
may originate from the nnmelted crnsts of partly differentiated planetesimals. Moreover, 
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Thomas et al. (2014) have suggested that the Massalia family may have olivine-rich objects 


among its predominantly S-type members. Accordingly, it is plausible that the X-type (135) 
Hertha is associated with the S-type Herthal family. 


In the following subsections, we consider four scenarios for the origin of the Herthal 
and Hertha2 families that might explain their different reflectance properties. None of the 
scenarios explain the observations completely; therefore the actual origin scenario may be a 
combination of these or other effects. 


5.1.1. Scenario 1: Unrelated families 


The possibility that both Herthal and Hertha2 originated from two unrelated collisions 
on different parent objects, one of which (Herthal) was catastrophically disrupted, cannot be 
ruled out here. One approach to understanding the type of collision expected to have formed 
the Herthal family involves comparison of the observed size-frequency distribution (SFD) of 
the family members with models of SFDs produced by catastrophic and cratering collision 
events. However, the incompleteness of the reflectance property databases, as well as the 
depletion of primarily smaller objects due to the nearby large resonances, have introduced 
signihcant bias into the family’s observed SFD. Future efforts will explore the connection 
between the Herthal family and (135) Hertha in detail, including numerical integrations to 
determine the possible evolution of the proper orbital elements of (135) Hertha (similar to the 


method used for (495) Eulalia by Walsh et al. 2013). In the meantime, the close dynamical 
proximity of both Herthal and Hertha2 families to the asteroid (135) Hertha remains strong 
evidence that they share this common parent, and warrants discussion of additional origin 
scenarios. 


5.1.2. Scenario 2: Single collision on a partially differentiated parent 

A single collision on a partially differentiated planetesimal could produce, in principle, 
a family of objects with varied material properties. The tighter Yarkovsky V of the Hertha2 
family {C = ±0.012 mAU vs. C = ±0.045 mAU for the Herthal family) would appear 
to suggest a younger age, and hence a separate collision. However, a single collision could 
still explain the two different Yarkovsky V’s if the tighter V reflects a decreased Yarkovsky 
semimajor axis drift rate instead of a younger age, with the objects of the Hertha2 family 
drifting outward more slowly after the original collision. 

A slower drift rate for the Hertha2 objects could result from their lower average albedo. 
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Py- The small numbers and small sizes of the Hertha2 family objects complicate the de¬ 
termination of the average albedo; however, it would have to be ~ 0.014 to account for 
the reduced semimajor axis spread, which is implausibly low (asteroid (135) Hertha has 
Py = 0.15). Thus, the relatively narrow V of the Hertha2 family cannot be accounted for by 
lower albedo. 


A higher average bulk density might alternatively explain the slow drift. The Yarkovsky 
drift rates scale inversely with density, so to account for the different V shape relative to the 
Herthal S-type family (with ps = 2.7 g/cm^), the average density that would be required 
for the Hertha2 objects would need to be 8.6 g/cm^. For comparison, the bulk density of 
Hertha is pn = 5.23 g/cm^, and the average for the Xk-type objects is pxy, = 3.8 g/cw? 
(Carry 2012). The larger density could be achieved if the family were composed primarily 
of metal objects; however, the family members’ apparent association with the lower-density 
(135) Hertha does not imply that this is the case. Thus, density differences probably cannot 
explain the reduced spread of the Hertha2 family, and a more recent formation event seems 
the most plausible explanation. 


Further evidence for a separate formation event comes from the observed ap — ep cor¬ 
relation among the Herthal family members’ orbits (Section]^, and the lack of one among 
the Hertha2 family members. If the Hertha2 members formed in the same collision as the 
Herthal members, they would be expected to share this correlation. Thus the most likely 
explanation is that two separate collision events occurred. 


5.1.3. Scenario 3: Two collisions on a partially differentiated parent 


A series of two collisions occurring on a partially differentiated parent could excavate 
families with varied material properties from different regions in the parent body. In the 
case of the Hertha families, an initial collision into a parent body’s unmelted, chondritic 
crust could form the S-type Herthal family, leaving a large X-type remnant, (135) Hertha. 
A subsequent impact onto (135) Hertha could then create the X-type Hertha2 family. 


A collision of magnitude sufficient to create the Herthal family might be expected to 
produce more than one remnant Hertha-like X-type among its fragments. The nearest object 
with reflectance properties similar to (135) Hertha is the large E-type (44) Nysa, which could 
have conceivably come from the original Herthal impact, although it is signihcantly displaced 
from the family in proper inclination. Zellner et ah (1977) and Kelley et ah (1994) attempted 


to link (44) Nysa and (135) Hertha in various differentiation structures, and Rivkin et ah 


(1995, 2000) demonstrated that the spectra of both objects share a hydration feature that 
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could suggest a common origin. Nevertheless, the issue of whether a hrst impact could have 
cleanly removed the S-type material leaving only one X-type remnant for subsequent impact 
remains uncertain. 


5.1.4- Scenario 4- Two collisions with shock darkening or space weathering effects 


The effects of shock darkening (through the presence of impact melt) can modify the 
apparent reflectance properties of asteroids without changing their basic material properties 
(Britt and Pieters 1994). In the case of the Hertha families, the impact that formed the 


Herthal family could produce remnants such as (135) Hertha with a signihcant amount of 
shock-darkened material, resulting in the later excavation of material with X-type appearance 
during the Hertha2 formation event. 


Reddy et ah (2014) demonstrated the effects of shock darkening among laboratory 


samples of LL chondrite material mixed with impact melt, and matched the resulting spectra 
with observations of the X-type Baptistina family members. In appearance, the difference in 
a* color between the Herthal and Hertha2 families (—0.07 vs. 0.13) is strikingly analogous to 
that of the nearby Flora and Baptistina populations: The Baptistina members lie in nearly 
the same place in a* and i — z color space as the Hertha2 members, as seen by comparison 


of the location of the Hertha2 family in Fig. 12 with the location of the Baptistina family in 


a similar plot (e.g.. Figs. 3e and 3f of Dykhuis et ah 2014); and both the Flora and Herthal 
populations show reflectance properties typical of S-type objects. 


Given the results of Reddy et ah (2014) for the Baptistinas, it is worth considering 


the possibility that the X-type appearance of (135) Hertha and the Hertha2 family members 
results from shock darkening, i.e. the presence of signihcant impact melt matrix surrounding 
fragments of the original S-type material. In this case, a plausible origin hypothesis is that 
the impact that formed the Herthal family left behind a shock-darkened fragment, (135) 
Hertha, which experienced a later collision that formed the darkened Hertha2 family. 

Another explanation for the difference in a* between the two Hertha families could be 
the effects of space weathering, which raises the a* value of exposed asteroid surfaces over 
time, resulting in higher a* for older families. In this case, the older Herthal family would 
be expected to show higher average a* values than the Hertha2 family, as observed. 

Consideration of the surface properties of the four subfamilies of the large, S-type Ko- 


ronis family (Molnar 2011) raises issues regarding both the shock-darkening and the space¬ 


weathering hypotheses. These subfamilies range in age from 5.8 My to 300 My, and display 


a linear trend in a* from 0.02 for the youngest family to 0.07 for the oldest (Molnar 2011 
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Dykhuis et al. 2014). This range of a* values is consistent with the color change expected 


as a result of space weathering, based on a study of various other families by Willman et ah 


(2010). The difference in a* between the two Hertha families is four times greater, so space 
weathering is unlikely to be the key factor in explaining this difference. If, instead, shock 
darkening is the dominant factor in the different appearance of the Hertha families, it raises 
the question of why similar effects have not led to more diversity among the Koronis families. 
The answer may simply lie in the unknown specific details of the collisions responsible for 
each family group. 


5.2. Eulalial and Eulalia2 


The low-albedo families in the Nysa-Polana region are thought to be the most likely 
source families for asteroid (101955) Bennu, the target of the OSIRIS-REx sample return 
mission (Campins et al. 2010| Walsh et al. 2013, Bottke et al. 2014). Bottke et al. (2014) 
determined the probability that Bennu originated from what we call the Polana family to 
be 70%, with the probability of an origin from what we call the Eulalial family of 30%. The 
Eulalia2 objects identihed here have similar orbital and reflectance properties, and might 
also be a candidate source population for Bennu. The possibility of a collisional connection 
between the Eulalial and Eulalia2 clusters, discussed below, could help further constrain the 
dynamics of the Bennu source region. 

The Eulalial and Eulaha2 clusters show very similar structure in ap — 77 space, and 
can be roughly fit with C envelopes that match, except for a semimajor axis offset of about 
6ap ~ 0.035 AU (Section 4.2, Fig. 15), with Eulalial centered at 2.46 AU and Eulaha2 


centered at 2.495 AU. These values bracket the current ap value for (495) Eulalia (ap = 
2.4867 AU). The Eulalial and Eulalia2 clusters also bracket (495) Eulalia in ep and sin ip 


(Fig. 14). Walsh et al. (2013) noted that the long-term evolution of (495) Eulalia ranges 


between ep values of 0.11 to 0.15; this range is contained within the larger range of the 
combined Eulalial and Eulalia2 families. 

In light of these dynamical coincidences (matching C envelopes, bracketing in each 
orbital dimension), it is plausible that the Eulalial and Eulalia2 clusters were originally part 
of the same collisional family, centered near the current position of (495) Eulalia. In this case, 
the twin “lobes” of the family protruding in opposite directions in ap,ep and sin ip could 
be explained by a non-isotropic ejecta pattern. In this scenario, due to the parent object’s 
proximity to the Jupiter 3:1 resonance, a signihcant number of the original objects would 
have been deposited directly into the resonance (though not beyond it, see the discussion in 


Section 2.5 of Walsh et al. (2013)). The resonance would then have dramatically reduced the 
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number of objects in the high-op, low-ep, low-sin ip lobe of the family. Subsequent Yarkovsky 
drift would then have removed another half of the remaining family (the prograde objects) 
via drift into the Jupiter 3:1. 


The above scenario could be tested dynamically, by modeling the best fit values of the C 
boundaries and op^o for fhe Eulalial and Eulalia2 clusters to confirm the similarity of C for 
the two and determine their offset in semimajor axis. This modeling is necessary to confirm 
the signihcance of the observed offset, as the sample of Eulalia2 members (from the magenta 
box in Fig. 14) is fairly diffuse and thus might be expected to contain a signihcant number 
of background objects. In addition, orbital integrators that incorporate the Yarkovsky and 
YORP effects could be used to simulate and reproduce the observed ap — H distributions of 
the family from a single collision. 


Lastly, spectral comparison of the reflectance properties of the Eulalial and Eulaha2 
structures can also test the above origin scenario. The SDSS and WISE databases provide 
reflectance data for a large number of asteroids in the region; however, the completeness 
of these databases is only around 20% and 10%, respectively. Future surveys such as the 
Gaia space mission of the European Space Agencjj^and the Large Synoptic Survey Telescope 
(LSST)[^are expected to signihcantly enhance these databases in the near future, providing 
opportunities for further analysis. The Polana and Eulalial and Eulalia2 structures are all 
similar in the wavelength range available to the SDSS, and a survey of 72 objects in the 
region shows spectral homogeneity in the NIR wavelengths as well (N. Pinilla-Alonso, in 
preparation). Future work in this area might possibly distinguish or further connect the 
Eulalial and Eulaha2 structures. 
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Fig. 1. —: (a) Density distribution in proper ep and sin ip of all objects with well-determined orbits 
between 2.1 and 2.5 AU. The red and blue boxes indicate the rough boundaries of the regions dominated by 
S- and C-type objects. The large asteroids (44) Nysa, (135) Hertha and (142) Polana are indicated by the 
red diamond, red square and blue diamond, respectively. The structure at ep ~ 0.16 and sin ip ~ 0.025 is 
the Massalia family; the Flora family muddies the upper left half of the plot. Subplots (b) and (c) show the 
distribution in proper semimajor axis and absolute magnitude of the objects in the blue and red boxes of (a). 
The gray and blue dotted lines mark the location of the Mars 1:2 and Jupiter 3:1 mean-motion resonances 
(ap = 2.419 AU and 2.495 AU, respectively). Several characteristic Yarkovsky “V” patterns appear in plots 
(b) and (c), suggesting the presence of multiple evolved families. 
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Fig. 2. —; Distribution of SDSS colors {a* and i — z) for all bodies included in Fig. for which adequate 
color data are available. The two major clusters seen here on either side of the boundary at a* = —0.015 
(dotted line) are also distinct in the three-dimensional parameter space {a* ,i — z,py). Generally, objects with 
larger a* have reflectance properties consistent with S spectral types, while those with smaller a* are darker 
spectral types (e.g., C- or X-types). 



Fig. 3. —: Proper orbital elements for all objects in Fig. with a* > —0.015. The large structure that 
occupies much of the box is here designated the Herthal family. As in Fig. the red diamond and square 
represent (44) Nysa and (135) Hertha, respectively. 

















27 



Fig. 4. —: Proper semimajor axis distribution for all of the objects within the red box in Fig. The 
dotted black, blue and gray lines indicate the locations of the Mars-crossing line and the Jupiter 3:1 and 
Mars 1:2 mean-motion resonances, respectively. The Herthal objects show a trend of decreasing eccentricity 
for higher semimajor axes. At the lowest semimajor axes and highest eccentricities, the distribution has been 
truncated by the Mars-crossing line; at the highest semimajor axes and lowest eccentricities, the distribution 
has been truncated by the Jupiter 3:1. The most pristine “boundary” of the family occurs for the objects 
with 2.3 AU < ap < 2.4 AU and ep < 0.185. The color gradient in eccentricity is used in Fig.lsl 
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Fig. 5.—: The population from Fig. plotted in op — i/ space, with color indicating eccentricity. Plots 
(a) and (b) are identical except for the ordering of the objects plotted: In (a), higher ep objects are plotted 
on top of lower ep objects; in (b), the lower ep objects are plotted on top to show the behavior at the lower 
eccentricities. The Yarkovsky “V” shape is evident among the family members. The black lines represent 
\C\ = 0.04 mAU, 0.06 mAU and 0.08 mAU; none of the lines ht the outer edge of the distribution well, due 
to the ep dependence which skews the distribution in ap. 
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Fig. 6.—: Equivelocity ellipses (i.e. the distribution in ep, sinzp of particles ejected isotropically at 285 m/s 
from the barycenter of the population) for various values of the true anomaly, /, of the collision event, plotted 
over the distribution from Fig. The best linear fit to the boundary, with slope m = —0.50 AU“^, is also 
shown (dashed line on left). 
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Fig. 7. —: The ep, ap distribution of the objects of the Herthal family, with the corrected values for proper 
semimajor axis (a^ calculated from Eq. [fusing slope m = —0.50 AU“^ and ep^o = 0.177), removing the 
skew in the distribution due to the initial collision. The boundary along the left edge is now vertical. The 
Mars-crossing region and Jupiter 3:1 resonance have noticeably depleted the objects at high ep (blue dots) 
and high ap (right side of the family), respectively. 
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Fig. 8. —: The C distribution of the Herthal family, computed from Eq. using the corrected Op with 
slope m = —0.50 AU“^. The curves represent the kernel density estimate (KDE) for the distribution, along 
with its first and second derivatives. The best-fit lower C boundary occurs at a local maximum in the second 
derivative, C = —0.045 mAU. 
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Fig. 9. —: The distribution of Herthal objects in a'p — H space. Objects with ep > 0.185 are not included 
to avoid the effect of the Mars-crossing region on the left edge boundary. The black lines shows the amount 
of semimajor axis drift that correspond to C{ap) values of ±0.045 mAU, with dotted lines representing the 
uncertainty of ± 0.005 mAU. The curves are fit to the left edge of the population; the curves on the right 
side are mirror images of those on the left for positive C values. 
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Fig. 10. —: Proper orbital elements for all objects in Fig. 0with a* <—0.015. Two structures appear within 
the blue box: the diffuse Polana family surrounding (142) Polana and the larger, denser Eulalial family near 
ep = 0.145 and sin ip = 0.047. Asteroid (495) Eulalia itself (green square) is offset in eccentricity from its 
nominal family. A cluster of low-a* objects (cyan box) is revealed near (135) Hertha, and is here designated 
“Hertha2.” The high-a* Herthal family totally obscures this cluster in Fig. [l] The faint structure at about 
sini = 0.025 is a low-a* remnant of the Massalia family, the spread of whose distribution in a* overlaps the 
a* = —0.015 boundary. 
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Fig. 11. —: The Hertha2 objects within the cyan box in Fig. 10 are plotted in op — space. They display 


a Yarkovsky “V” bounded by C = ±0.012 mAU (cyan lines). The C parameter boundaries for the Herthal 
family from Fig. are plotted as black lines for reference. Asteroid (135) Hertha is a plausible parent of the 
Hertha2 family as well as of Herthal, although both families have distinct C envelopes as well as reflectance 
properties (Fig. 12). 



a* color 


Fig. 12. 


The same as 


Fig-i 


but restricted to the range of ep and sini given by the cyan box in Fig. 


10 


The Hertha2 family (near a* = —0.07) is distinct from both the larger Herthal family (near a* = 0.13) and 
the Polana/Eulalia families (near a* = —0.13). 
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The low-a* objects within the orbital range plotted in Fig. 10 show at least two Yarkovsky 


Fig. 13. 

V signatures, as noted by Walsh et al. (20131 (cf. their Fig. 2). One structure is roughly centered on (142) 
Polana (blue diamond); the other is roughly centered on (495) Eulalia (green square). The curves correspond 
to C = ±0.169 mAU (blue) and C = ±0.092 mAU (green); the cyan lines show the Yarkovsky envelope for 
the Hertha2 family (cf. Fig. 11), which is included in this low-a* population. 



Fig. 14.—: Orbital distribution of low a* objects, similar to Fig. 10 but with adjusted ranges in ep and 


sin ip to include the region near (495) Eulalia (green square). A diffuse cluster appears within the magenta 
box. Walsh et al. (2013) associated this grouping with the object (112) Iphigenia (green diamond). The 
green box indicates the core of the Eulalial family, which makes up much of the population above the green 


line in Eig. 13 
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(b) 


Fig. 15. —: (a): The same as Fig. 13 but with objects with orbits in the green box in Fig. 14 (predominantly 


Eulalial members) highlighted in green. These objects are roughly centered on (495) Eulalia and bounded 
by the model curve of C = —0.092 mAU. (b): The same as (a), but with overplotted objects from within the 
magenta box in Fig. 14 These objects are roughly centered on (495) Eulalia as well, but do not extend all 
the way to the green curve, thus indicating that the structure within the magenta box in Fig. is distinct 
from the Eulalial family, and hence designated Eulalia2. 








